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ABSTRACT 

Based on a new estimation of their thickness, the global properties of relativistic slim accretion disks are in- 
vestigated in this work. The resulting emergent spectra are calculated using the relativistic ray-tracing method, 
in which we neglect the self-irradiation of the accretion disk. The angular dependence of the disk luminosity, 
the effects of the heat advection and the disk thickness on the estimation of the black hole spin are discussed. 
Compare to the previous works, our improvements are that we use the self-consistent disk equations and we 
consider the disk self-shadowing effect. We find that at the moderate accretion rate, the radiation trapped in 
the outer region of the accretion disks will escape in the inner region of the accretion disk and contribute to 
the emergent spectra. At the high accretion rate, for the large inclination and large black hole spin, both the 
disk thickness and the heat advection have significant influence on the emergent spectra. Consequently, these 
effects will influence the measurement of the black hole spin based on the spectra fitting and influence the 
angular dependence of the luminosity. For the disks around Kerr black holes with a — 0.98, if the disk inclina- 
tion is greater than 60°, and their luminosity is beyond 0.2 Eddington luminosity, the spectral model which is 
based on the relativistic standard accretion disk is no longer applicable for the spectra fitting. We also confirm 
that the effect of the self-shadowing is significantly enhanced by the light-bending, which implies that the non- 
relativistic treatment of the self-shadowing is inaccurate. According to our results, the observed luminosity 
dependence of the measured spin suggests that the disk self-shadowing significantly shapes the spectra of GRS 
1915+105, which might lead to the underestimation of the black hole spin for the high luminosity states. 

Subject headings: accretion, accretion disks - binaries: close - black hole physics - X-rays: stars 



1. INTRODUCTION 

The standard acc retion disk model (SSD, 
IShakura & Sunvaevl 119731) an d its relativistic generaliza- 
tion dNovikov & Thorndl 19731) are widely used in modeling 
the spectral energy distribution (SED) of both active galactic 
nuclei (AGNs) and the galactic black hole binaries. Recently, 
based on the relativistic SSD, the spins of the black holes 
in several galactic black hole binaries are measur ed by 
fittings of their SEP dLi et all l200l Shafee et al. 2006; 
iMcClintock et al.ll2006t iMiddleton et alJ 120061) . One major 
advantage of standard accretion disk model is its simplicity, 
i.e., it assumes that the fluid undergoes the Keplerian motion 
and the energy generated by the viscous heating is radiated 
locally. Besides, in this model, the disk terminates at the 
Innermost Stable Circular Orbit (ISCO), and the matter inside 
ISCO is cold. 

However, when the accretion rate is relatively high (M > 
O.lMcddX m the inner region of the accretion disks, several 
effects will make these assumptions invalid. One such effect 
is the heat advection. When the vertical diffusion time scale is 
larger than the accretion time scale, substantial amount of heat 
is trapped i nside the accretion flow and not able to esc ape to 
the infinity jBegelmanlll978l:lBegelman & Meierj|1982l) . This 
will make the local flux different from that from the SSD 
model and affect the profile of the radial temperature. Also, 
the advected heat will provide an additional pressure, which 
makes the disk rotation no longer Keplerian. Besides, because 
the heat trapped inside the accretion flow, the region inside 
ISCO is not cold any more. 
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The self-consistent treatment of these effects needs to 
model the energy and momentum balance in the disk, one 
widely used method is based on the radiative MHD simu- 
lation. Although the numerical simulations have the advan- 
tage of treating the physical effects self-consistently, the huge 
amount of computation makes it difficult to give the quali- 
tative fitting of the observational data. One complementary 
approach is to solve the energy and momentum equations of 
the accretion flow in terms of the ordinary differential equa- 
tions (ODEs). Compared to the simulations, this method 
can effectively take into account the dynamical effects, while 
maintaining its concis eness. In this work, we use the ODE 
form disk equations (iPaczvnski & Bisnovatvi-Koganl 11981b 
Abramowi cz et all 119881: IChen & Taaml 119931: iPeitz & Appll 
1997; Gammie & Pophaml ll998l) . and explore the solutions 
at both the medium a nd high accretion rate ( where the disk is 
call the "slim disk"; Abra mowicz et al|1988l). We use the rel- 
ativistic disk equat i ons in our work (lAbramowicz et alii 19961 : 
Beloborodov 1998; Shimura & Manmoto 2003) to investigate 
the effects of the black hole spin. 

The physical properties of the black hole and the accre- 
tion disks are deduced mainly from their observed SED. It is 
well known that the structure and the SED of the slim disk 
differ significantl y from that of SSD (IWang & Zhod [19991: 
IWang et alJfl999l) . The disk SED is influenced by a handful 
of effects, such as the Doppler effect, the relativistic beam- 
ing, the gravitational redshift, the gravitational -bending, and 
the disk self-shadowing. As the accretion rate increases, 
the disk b ecomes geometrica lly thick, and the disk self- 
shadowing (IWu & Wangl 2007) becomes signific ant. In this 
work , we employ a ray -tracing approach (e.g. iFanton et al] 
ll997l;ICadez et al.ll 19981) to explore the final spectra at the dif- 
ferent viewing angles. 
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Due to their disk-like geometry, the radiation from the ac- 
cretion disks is far from the isotropy, and the emergent lumi - 
nosity is thus angular dependent (e.g. Hein zeller et al.l l2006). 
For a slab of optically thick material without velocity, the disk 
luminosity simply varies as L ~ cos 6 in the fiat spacetime. In 
fact, because of its thickness and rotation, the angular depen- 
dence of the disk luminosity is expected to be different from 
that of the simple slab, further analysis is thus needed. In this 
paper, we calculate the angular dependence of the disk lumi- 
nosity under different inclinations and the black hole spin. 

The spin of GRS 1915+105 is debated in th e literature 
(McClintocketaLl 120061: iMiddleton etaLl 120061). B y fit- 



ting the disk blackbody component, McClintock et al. (2006) 
found that the spin is near extreme, while Middleto n et al.l 
(2006) found the spin is moderate, with a ~ 0.8. To show 
how the spin determination is influenced by the heat advec- 
tion and disk s elf-shadowing, we re-sim ulate the observation 
carried out by iMcClintock et al.l ([2006) with our model and 
confirm the luminosity dependen ce of the measured spin as 
found by IMcClintock et al.l (120061) . 

In this work, the black hole mass in our canonical model 
is taken to be 10M Q . However, due to the physical similari- 
ties between the accretion disks around the stellar-mass black 
holes and those around the supermassive black holes, our re- 
sults can be easily generalized to the supermassive black holes 
and predict the SED of AGNs under the simple scaling. 

This paper organized as follows. In section 2, we de- 
scribe briefl y our d isk m odel and the ray-tracing method. In 
Sectior inT Section 13.31 we present our numerical results of 
the global disk structure. In Section O we present our results 
on the SED of the disk, focusing on the differences between 
the SED of the dynamical disk model and that of the SSD. The 
angular dependence of the total luminosity is given in Section 
[5] In Section [6] implications for estimating the black hole 
spin are discussed. In Section [7] various observational impli- 
cations of the disk self-shadowing as well as the spin of GRS 
1915+105 are discussed. Finally, in Section [8] we summarize 
our conclusions. 

2. THE RELATIVISTIC SLIM MODEL 

2.1. The Dynamical Disk Model 

The basic equations governing the disk structure used in 
our wo rk are similar to those given in Shi mura & Manmotol 
(2003). In the Boyer-Lindquist coordinates, expanding 
around the equ atorial plane up to (z/r) ° terms, the metric 
takes the form ( Novikov & Thornell 19731 the geometric unit 
G = C = M = 1 is taken) : 

ds 2 = ~r 2 ^-dt 2 + 4(# - udt) 2 + ^-dr 2 +dz 2 , (1) 



with 



A = r - 2r + a , 



A = r i + r 2 a 2 
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where a is the dimensionless spin of the black hole. 
The mass conservation takes form 



M = 27rA 1/2 E 7r t/ , 



(2) 



where M is the mass accretion rate, and W and £ are inte- 
grated pressure and density, defined as W = f p dz, £ = 



/ p dz. Here, p means pressure and p means density, the sub- 
script "0" denotes the corresponding value at the equatorial 
plane, and H is the scale height of the accretion disk. j r is 
the radial Lorentz factor of the fluid. 
The angular momentum conservation reads 



M (T T ^ 
Ztt 



A 3/2 A l/2 3 



-W, 
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where L is the specific angular momentum, L- ln is the spe- 
cific angular momentum at the inner boundary of the accretion 
fl ow, and 7^ i s the L orentz factor in the <f> direction as defined 
in Manmoto (2000). As the vi scosity prescription used here 
is conventional (Naravan 1992), at the inner part of the accre- 
tion disk, our result may be inaccurate in principle. However, 
the causal viscosity is still somewhat putative and differs from 
the conventional one only close to or inside the ISCO. Thus, 
the disk spectrum in our Kerr case is robust in the sense that 
the sonic point is very close to the event horizon of the black 
hole, while the disk spectrum in the Schwarzschild case may 
be shifted a bit due to the different treatment of the prescrip- 
tion of the viscosity. 

The momentum equation and energy equation take the 
forms as follows: 
Momentum equation: 



2 V dV 



i aw _ i\a (n-n+){n- 

E~dr~ ~ r^A fi+f^ 



(4) 



where 
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i 



r 2 ± a 



is the Keplerian angular velocity of the direct ("+ sign") and 
retrograde ("-" sign) circular orbit, respectively. And here 

r 2 L 2 \V2 



70= (1 
is Lorentz factor of v$, 



= uj + Q = uj - 



Energy equation: 



Qvis ~ Q&dv + 



ad ' 



(5) 



here, Q^ is is the viscous heating rate per unit surface, F d is 
the radiative cooling term, and Q~ dv is the effective cooling 
induced by heat advection. They can be specified as follows: 



Q + 

^ VIS 



47rr A 1 / 2 r ^ m ^ dr 
M W 1 sdlnW 



* adv 



47rr E T3 — 1 V dr 

dlnS ^ d\nH 
1 1— (1 l + 1J— -, 

dr dr 

where Ti and T3 are defined in lKato et alj (fl998l) . and 

_ 8acT 4 
rad ~ 3np H ' 
The Rosseland mean of total opacity can be expressed as 

K = 0.4 + 0.64 x 10 23 pT~ 7/2 . 



(6) 



(7) 



(8) 
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The total pressure can be expressed as the sum of the gas 
and radiati on pressure. For simp l icity, w e drop the numerical 
factors in Shim ura & Manmotol ((2003), and the equation of 
state is given by: 



(9) 



where T is temperature of the disk. Another equation is the 
vertical thickness equation, which reads 



H 



(10) 



This is a mod i ficatio n to the equations of 
Shimura & Manmot o (2 003), and its derivation will be 
presented in Section l2T2l 

The above equations are transformed into a set of first-order 
differential equations with variables v T and W, and then inte- 
grated from the outer boun dary at about 10 5 r g . 

The transonic solutions (Liang & Thompson 1980) are ob- 
tained using a shooting method, in which we adjust L ln re- 
cursively until the smooth transonic solution is found. During 
the solving process, despite that the factor g pa takes form 



•9i 



pa 



7^ 



(r 2 + a 2 ) 2 + 2Aa 2 



Aa 2 



(r 2 + a 2 ) 2 



(11) 



a simplified version of factor g pa is used when finding the 
global transonic solution: 
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(r 2 
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However, the location of the photosphere (Peitz & App lll997h 
is still determined according to Equation ([TJl l. When this sim- 
plification is taken, at the inner regions of the accretion disks, 
the resulting disk photosphere and disk thickness no longer 
agree with each other. But the difference is not significant. 
We thus neglect this inconsistency and use the location of the 
disk photosphere in the ray-tracing calculation. 

2.2. Estimation of the Vertical Thickness 

The vertical thickness is a major uncertainty in the accretion 
disk theory, and is claimed to be dependent on various phys- 
ical assumptions such as the vertical dissipation profiles and 
the interplay between the relaxation time scale and the free- 
fall time scale. However, Sadowski et al] (120091) showed that 
the vertical thickness of an accretion disk in the non-advective 
regime can be approximated using this simple formula: 



kuT, 



(13) 



where k is the total opacity, a is Stefan-Boltzmann constant 
and z is the photosphere location. T e g is the effective tem- 
perature defined as T e g — (F^/cr) 1 ^ 4 and = r~ 3 / 2 . 
Equation (fT3l simply means that at the disk photosphere, the 
vertical radiation pressure balances the gravity of the central 
black hole. In this work, despite the fact that Equation ( fT3l 
has only been shown to be valid in the non-advective regime, 
we also use it in the advective regime, since the physical ar- 
gument is still valid. With this assumption, we can get the 
expression for the location of the photosphere: 



3 ^rad ri _2 
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Another constraint comes from the vertical radiation flux. 
According to the diffusion approximation, assuming the dom- 
inance of the radiation pressure, the radiation flux can be ap- 
proximated as 



-^rad — 



8acT 4 
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(15) 



We thus make the physical requirement that the thickness z 
obtained in Equation < TT~4T > should be consistent with the thick- 
ness H in Equation (T5[ . Combining these two equations, we 
obtain the expression for modified disk thickness: 



H 



(16) 



It should be noted that Equation ( [T6l l is accurate when the disk 
is radiation pressure dominated. At the outer regions of the 
accretion disks, where the gas pressure dominates, Equation 
(16) deviates from the previous one by a factor of y/12 (the 

I/O 

previou s expression takes form H = Cs/OkS'pa ). How- 
ever, as [Sadowski et al. (2009) showed, the thickness in the 
outer region of the accretion disks is not well constrained, 
and it is dependent on the dissipation profile in the verti- 
cal direction. As the contribution of radiation from the gas- 
dominated region is relatively unimportant and the final emer- 
gent SED is unaffected by our modifications significantly, we 
use EquationdToT) throughout the paper. One concern about 
the estimation of thickness is that t he effect of magn etic field 
may influence the disk thickness. iHirose et alj (120061) pointed 
out that the location of photosphere of the accretion disk is 
modified by the magnetic field significantly. However, in 
our case, the shadowing effect happens at the inner region, 
where the radiation pressures dominates. Thus the thickness 
and the self-shadowing might not be significantly affected by 
the magnetic effect. The magnetic field might also change 
the inner boundary cond ition significantly (e.g. Pen na et al] 
120 10t iNoble et"aTll2010l and references therein). To inves- 
tigate the effects of the magnetic fiel d, we should do the 
general relativistic MH D simulations (Gammie et al.1 120031 : 
iDe Villiers et alj [2003h . which is beyond the scope of this 
work. 

2.3. Ray- Tracing Method 

After obtaining the global structure of the relativistic disk, 
we use the ray-tracing technique to calculate the emergent 
spectra with inclusion of all the relativistic effects. The 
ray-tracing metho d is based on the analytical integration of 
the null-geodesics .(IChandrasekharl 19831: iRauch & Blandford 
Il994t ICadez et alJI 19981) , and th e subroutines are the same as 
those used in I Yuan et al.1 (120091) . 

One major improvement compared to earlier work (e.g. 
Li et al. 2005) is that we include the disk thickness in the ray- 
tracing calculation. To achieve this, we start from the infinity, 
and then search along the each photon trajectory for the inter- 
section point. 

Once the intersection points are found, we can calculate the 
local temperature of the emitting points in our disk model and 
calculate the redshift factor g. The emergent spectra are then 
obtained by integrating over the photometric plane according 
to 

F Vo =fflJ^, d7) 



(14) where a and /3 are two impact parameters, and the redshift 



factor is 



9 = — 



(18) 



here D is the distance of the object expressed in geomet- 
rical units and I v takes t he diluted blackbody form (e.g., 
Shi mura & Takaharalfl995h . 
The redshift factors are calculated as follow: 



^obs 



server lr=oo, 0=0 



P^ U F\md\ r = 



(19) 



where p u is the four-momentum of the photon, given by 



wnere p u is 
ICarterl (1968) 



(20) 



The four-velocity of the fluid is given bv lYuan et all (120091) 

u Fiuid = (lrl^e~",'y r (3 I e~' Jl1 , , 7 r 7 (a;e _l/ + /V"^)) , 

(21) 

where e v = (SA/A) 1 / 2 , = {A sin 2 djY,) 1 ! 2 . Thus we 
have the expression for redshift factor: 
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A 1 / 2 
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r 2 A l/2 



(22) 



where "em" means the emitting points, and the "+" sign for 
the trajectories of the photons without r-turning point and "— " 
sign for those with r-turning point. In Equation(20), 

R = r 4 + {a 2 -\ 2 -Q)r 2 + 2(Q + {\-a) 2 )r-a 2 Q . (23) 

here the definition of A and Q are the same as that of 
I Cunn ingham & Bardeenl (119731) . 

In the ray-tracing calculations, the outer boundary condi- 
tion of the integration is set at 10 5 r g . In the cases when 
the emitting radius is so large that the disk solutions do not 
cover, we use the results from SSD instead, because the differ- 
ences between the transonic disk model and SSD are indeed 
insignificant at the large radii. 

One limi tation in our calculation is that the returning ir- 
radiation (Ivan Paradiis & McClintockll 19951: ID" et al.ll2005l) is 
neglected. In the inner regions of the accretion disk, due to 
the strong gravity, the radiation emitted at one part of the 
accretion disk may reach other parts, which may result in 
a modification to the disk structure and the emergent spec- 
tra. Certain disk geometry may e nhance such kind of effect 
(van Paradiis & McClintock 1995). In the case of the thin ac- 
cretion disk, this effect has moder ate influence ( about several 
percents) on the emergent spectra dLi et alj2005l) . In our case, 
due to the thick geometry, since the disk is shadowed by itself, 
the self-irradiation effect may be stronger. As the inclusion of 
that effect requires significant amount of calculation, for sim- 
plicity, we neglect it in our current work. 

Another limitation in our calculation is t hat the effect of 
limb-darkening has not been considered. As iLi et ail (2005) 
have pointed out, the limb darkening has effect on the emer- 
gent spectra, especially when the inclination angle is large, 
and the deviation induced by this effect is ab out 10% at low 
energy band when the disk is viewed at 85° dLi et alj|2005l) . 
In this work, this effect has not been included. According to 
the limb-darkening law Ig ~ 1/2 + 3/4 cos 6, our results are 
still robust at relatively small inclination angles. 



3. STRUCTURE OF THE RELATIVISTIC DISK 

We present our results on the structure of the accretion disk. 
Throughout this paper, the critical accretion rate is taken as 



M edd = 16- 



J cdd 



(24) 



where the numerical factor comes from the energy generation 
efficiency of the Schwarzschild black hole. For the sake of 
comparison, we adopt the same critical accretion rate for both 
the Schwarzschild and Kerr black holes. This is d ifferent from 
the convention used in iMcClintock et al.l (120061). Similar re- 
sults have already been obtained by Sadowski (2009), and one 
major difference between our work and theirs is that our work 
is based on a different estimation of the vertical thickness. 

In our canonical model, the mass of the black hole is taken 
to be 10A/q, while the viscosity coefficient a is set to be 0.1. 
The dimensionless accretion rates are taken to be 1 /4, 1 /2, 1 
and 2 respectively, ranging from the sub-Eddington systems 
to the moderate super-Eddington ones. A special case for the 
low accretion rate M = 1/16 is also discussed when a = 
0.98. 

3.1. Disk Dynamics 

Figures Q][3] show the dynamical properties of our accretion 
disks. The upper panels in Figure[T]show the radial velocity as 
a function of radii, with the different accretion rates. Because 
we use fully relativistic equations, near to the event horizon 
of the black hole, as expected, the radial velocity is very close 
to the speed of light. At the small radii, for the same accretion 
rate, the higher the spin, the smaller the radial velocity. At 
the large radii, the radial velocity of the accretion flow is rel- 
atively small, and it is insensitive to the black hole spin. This 
is because that at the large radii, the disk structure resembles 
that of SSD. With the increase of the accretion rate, the radial 
velocity becomes larger. 

The lower panels of FigureQ]show the radial distribution of 
the surface density under the different accretion rates. Gener- 
ally speaking, at the small radii, the surface density is small. 
These results can be simply understood in terms of the non- 
relativistic form of the matter conservation: 

M = 2irrEv t , 

where E is the surface density of the accretion disks, and v r 
is the radial velocity (the relativistic one is similar but with 
an additional corrections related to the radial Lorentz fac- 
tor 7 r and the black hole spin a). Since the accretion disks 
around Schwarzschild black holes tend to have higher radial 
velocity, their surface density at given radius is smaller. An- 
other noticeable feature is that the density hop occurs when 
M = 0.25M c dd, a = 0, which is accompanied by the veloc- 
ity decrease and re-increa se in the upper panels. This pe culiar 
behavior is not shown in Shimura & Manmoto (2003) since 
they did not explore the disk solution at such a low accre- 
tion rate, but can be found in the p revious papers solvin g the 
non-relativistic slim dis k structure dChen & Wang 2004|) and 
a recent relativistic one (Sadowski 2009, with explanation). 

Upper panels of Figure |2] show the rotation velocity of the 
accretion flow. The rotation velocity is insensitive to the ac- 
cretion rate, especially at the outer region. To show how the 
fluid motion deviates from the Keplerian motion, the ratio of 
the rotational velocity to the rotational velocity of the Keple- 
rian motion / f? cp are shown in the lower panels of Figure 
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a=0 a=0.98 




r(r g ) r ( r 9 ) 



FIG. 1. — Radial velocity and surface density as a function of radii. From thin to thick, the lines stand for the solutions for the accretion rates l/16M c( jd, 
l/4Af cc id, l/2A;f c[ jd, lAif e< jd and 2M dd respectively. The left panels correspond to Schwarzschild black hole, and the right panels correspond to Kerr black 
hole with a = 0.98. The black hole mass is taken to be 10-Mq. 





Fig. 2. — Rotation velocity and the ratio of rotation velocity to Keplerian rotation velocity v^/v^^ vs. radii. The symbols of the lines are the same as in 
Figure[T] The left panels correspond to Schwarzschild black hole, and the right panels correspond to Kerr black hole with a = 0.98. 
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12 When the accretion rate is low, the rotation velocity ratio 
is very close to unity, which means the assumptions made in 
SSD are accurate. When the accretion rate is relatively high, 
the rotation of the fluid becomes sub-Keplerian. In some re- 
gions, the rotation is super-Keplerian. Although there exists 
some differences, the deviation is no more than 10% at several 
gravitational radii, which means the assumption of the Keple- 
rian motion is still va lid at the most of the regions of the a c- 
cretion disks (see e.g. Igumenshchev & Abramowicz 2000). 

Figuere[3]shows the effective temperature T e g as a function 
of radii. Here, T e ff is defined as T c ff = (-^'rad/ ') 1 ^ 4 ' where 
F r ~, is the radiation flux per unit area and a is the Stefan- 
Boltzmann constant. Similar to the above figures, the differ- 
ences between the Schwarzschild and the Kerr black holes are 
all very small at large radii. With the increase of the accre- 
tion rate, the temperature increases. When the accretion rate 
is high enough, the temperature is still high even inside the 
ISCO, which suggests that, contrary to the assumptions made 
in the SSD model, the matter inside the ISCO is not cold. 

3.2. Disk Thickness 

Figure|4]shows the thickness of the disk vs. the radii. In or- 
der to give a instinctive picture of the disk thickness, we use 
linear coordinate, and show only the results in the inner part 
of the accretion disk. This thickness profile is thus correspond 
to the surface of the accretion disk. The imprints of ISCO ap- 
pear in the disk thickness profile. This feature is particularly 
prominent at low accretion rate, wh ere the disk surface agrees 
wit h the simulation carried ou t by lHawlev & Krolikl (2002) 
and iRevnolds & Fabian] (120081) . When the accretion rate is 
high enough, the imprints of ISCO disappear, and the disk 
shows a continuous profile down to the event horizon of the 
black hole. 

However, in some cases, the estimation of thickness may be 
incorrect. Since the estimation of the thickness is based on the 
vertical hydrostatic equilibrium, it depends on the assump- 
tion: 

a F J uid < 9z , (25) 

which means the acceleration in the vertical direction could 
never exceed the acceleration provided by the vertical gravi- 
tational force (or that the disk has time to relax to the hydro- 
static equilibrium). We checked our solutions using this for- 
mula, and found that in some cases such as M = 0.25, a = 
Equation d25l ) is not satisfied. Thus the disk thickness may be 
under-estimated at the regions near ISCO, in some cases. 

3.3. Energy Dissipation Inside ISCO 

In Figure [5] we present our results on the global energy 
balance of the accretion disks. The energy conservation in the 
accretion disks takes the simple form as follows: 

Qt is = Km + Qadv , (26) 

where Q+ is is the viscous energy generation rate, F r ~ d is the 
cooling term due to radiation from the accretion disk, and 
Q~ dv is the effective cooling induced by the heat advection. If 
Q^ is ~ FZ<i> the disk is very similar to the standard accretion 
disk since the energy generated by the viscous processes is 
radiated locally. If Q^ is ~ Q^dv' most °f me heat generated 
is trapped inside the accretion disk. 

Figure [5] shows the ratio of the emitted radiation from the 
disk surface to the total heat generated by the viscous pro- 
cesses. At about 20 r„, the heat advection works, as a result, 



the ratio is less than 1, However, in the inner region, in many 
cases, F~ d dominates over Q^ is , which means that substan- 
tial amount of energy is still being radiated away. The source 
of the energy mainly comes from the energy trapped in the 
region outside the ISCO. Thus our results suggest that the ra- 
diation trapped outside the ISCO can be allowed to escape at 
the regions close to the black hole. This kind of mechanism 
is new, since different from the advective cooling, in some 
cases, the effect of advection will heat the disk up, resulting 
in "advective heating". 

Besides from getting self-consistent disk solution, we find 
that when the accretion rate is low enough such as M ~ 1 /4 
for a = and M ~ 1/16 for a = 0.98, the SSD model is 
valid. 

4. EMERGENT SPECTRA 

We present the resulting emergent spectra based on the ray- 
tracing method. We assume that the local emission of the 
disk takes the diluted b lackbody form w i th the spectra harden- 
ing factor e quals 1.7 dRoss et al.lll992t IShimura & Tak ahara 
1 1 9931 119951). which is actually dep endent on disk l umino sity 
as shown in Davis et al. (2005) and lDone & Davis (2008b. 

To show the effect of the photon trapping and the disk thick- 
ness, we consider three cases. In the first case (case A), we 
calculate the spectra from the relativistic standard accretion 
disks. In this case, our spectra are completely identical to the 
spectra generated by KERRBB with exactly the same parame- 
ters. In the second case (case B), we calculate the spectra from 
the relativistic slim accretion disks, assuming that the disk is 
still geometrically thin. This is for the purpose of showing 
the effect of the heat advection, the sub-Keplerian motion, the 
radial plunging and the radiation from the plunging region. In 
the third case (case C), based on case B, we consider the effect 
of disk thickness to show the disk self-shadowing. 

The emergent spectra are shown in Figure [6] The differ- 
ences between case A (black solid line) and case B (blue 
dotted line) are due to the combination of the effects of the 
photon trapping, the radiation from the plunging region and 
the motion of the fluid. The effect of the motion of the fluid 
on the emergent spectra is determined mainly by the Lorentz 
factors, 7 r , 7^. At the <fi direction, the velocity differences be- 
tween the slim disk model and the thin, Keplerian SSD model 
is relatively small, thus the effect is not significant. At the 
r direction, the radial velocity works. Since matter plunges 
into the black hole, in most cases, this will introduce a net 
redshift. At some rare ca ses where the photon trajectorie s 
have r turning points (e.g Cadez, Fanton, & Calvani 1998), 
the radial velocity field will introduce net blue-shift. But this 
blue shift is usually unimportant since this kind of trajectory 
is rare. The radial velocity effect is significant only near or 
inside the ISCO. 

The other important effect is the photon trapping. Since the 
previous works focus on the application of slim disk model in 
the highly accreting systems, the effects of the photon trap- 
ping in their cases are thus to bring down the luminosity, and 
highly attenuate the flux at the high energy band. However, 
for the relatively low accretion rate, the photon trapping be- 
haves quite differently. In fact, contrary to the previous view 
that photons are trapped and not allowed to radiate away, sub- 
stantial amount of trapped photons is allowed to re-radiate in 
the plunging region where the energy dissipation is relatively 
insignificant. This effect will in fact distort the spectrum and 
hence affect the spin measurement. This phenomenon is com- 
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FIG. 4. — Vertical thickness as function of radius. The symbol of the lines are the same as in Figure[T] The left panels correspond to the Schwarzschild black 
hole, and the right panels correspond to the Ken' black hole with a = 0.98. 




FIG. 5. as function of r. The symbol of the lines are the same as in Figure|4] 
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mon for both the non-rotating and rotating black hole systems 
when the accretion rate is moderate. 

The differences between case B and case C are due to the 
inclusion of the thickness of the disk in our calculations. It is 
clearly seen that the inclusion of the effect of the disk thick- 
ness will bring down the flux, especially at high energy band, 
and is more prominent at the higher accretion rate, higher in- 
clination angle and larger spin. When the accretion rate in- 
creases, the disk thickness increases, and the self-shadowing 
happens at a relatively small inclination angle. For the Kerr 
black hole, since the energy generation rate is higher, the disk 
self-shadowing effect is more prominent. 

Since the differences between case A and case B are what 
referred to as the photon trapping effect, the differences be- 
tween case B and case C are due to the effect of the disk 
self-shadowing, which has never been worked out before. In 
fact, when the accretion rate is high, the effect of the disk 
self-shadowing is more significant than that of the photon- 
trapping effect alone. Thus, in order to study the spectra from 
the super-critical accretion disks, we have to consider those 
two effects simultaneously. 

According to our results, for the disk around a Kerr black 
hole with a = 0.98, if the disk inclination is 60°, at about 0.23 
of Eddington luminosity, the effects of the heat advection and 
disk self-shadowing begin to be important. Although in our 
case (the bottommost set of curves in Figure [6j, the effects 
of the heat advection and disk self-shadowing tend to cancel 
each other, our results set a critical luminosity beyond which 
the relativistic standard accretion disk model is no longer ap- 
plicable for the spectra fitting. 

5. ANGULAR DEPENDENCE OF THE LUMINOSITY 

We explore the angular dependence of the disk luminosity. 
In the case of a thin slab of the stationary gas, the angular 
dependence of emergent luminosity should be L ~ cos 8. In 
the case of an accretion disk, the angular dependence of the 
disk luminosity L(6) is also influenced by several other ef- 
fects, including the disk geometry and the disk rotation. In 
Figures |7][8] we show both the luminosity in the X-ray band 
(0.1 ~ lOOKeV), which is of observational interest, and the 
total luminosity. 

Figure Q shows the luminosity in the X-ray band as a func- 
tion of the inclination. The upper panels show the results for a 
non-spinning black hole, and the lower panels show the spec- 
tra for a spinning black hole with a = 0.98. For the standard 
accretion disk model, the angular dependence of the luminos- 
ity deviates from L ~ cos 9 due to the Doppler shift induced 
by motion of the fluid. For the disk around a Schwarzschild 
black hole, the effect of the Doppler shift tends to flatten 
the profile, which means that the radiation is more weakly 
beamed than that of a single slab, for the disk around a Kerr 
black hole, this effect is so strong that the angle at which lu- 
minosity peaks is not 0°. 

The dotted lines in the subpanels of Figure|7]show the angu- 
lar dependence of the X-ray luminosity for a dynamical thin 
accretion disk which is different from the standard accretion 
disk since its velocity field is different and the photon trap- 
ping effect is taken into account. For the low accretion rate, 
as we have discussed, the effect of photon trapping is compli- 
cated. For the high-enough accretion rate, the effect of heat 
advection will highly attenuate the total luminosity. 

The solid lines in the subpanels of Figure [7] show the an- 
gular dependence of X-ray luminosity for a dynamical thick 
accretion disk. Compared to the blue dotted lines, the effect 



of disk thickness is taken into account. In all cases, the inclu- 
sion of disk thickness tends to attenuate the radiation, even in 
the face-on case, since light is bend by the strong gravity. An- 
other difference caused b y the disk thickne ss is a mild beam- 
ing of radiation (e.g. lHeinzeller et al.l2006l) . which means that 
the disk is more luminous when viewed face-on than viewed 
edge-on. 

The angular dependence of the total luminosity is shown in 
Figure[8] The trend is similar to the previous figure. However, 
the effect of the motion of the fluid is relatively weak. The 
reason is that in calculating the X-ray luminosity, the low- 
energy photons can be Doppler boosted into the X-ray band 
and contributed to the X-ray luminosity. However, that effect 
will not significantly influence the total luminosity. 

The mass of the black hole in our model is taken to be 
IOMq. Since the disk thickness in the inner regions of the 
accretion disks is dependen t on the accretion rate but inde - 
pendent of black hole mass dKato et al.lll998tlWatarail 2006). 
for a given 9 cl -n, the effect of self-shadowing is only a function 
of the black hole accretion rate. Given the same Eddington- 
scaled luminosity, since the other physical parameters such 
as T (the temperature or equivalently the characteristic fre- 
quency) and L to t (the normalization of the flux) scale as 

T~M~i (27) 

itot ~ M , (28) 

our results can be scaled to a considerable range of parame- 
ters, and then can be used to fit the observational results of the 
supermassive black holes. 

6. IMPLICATIONS FOR SPIN MEASUREMENT 



TABLE 1 
Parameters of KERRBB 



Parameters Name 


Description 


Value 


Freeze 


V 


Inner torque 





Yes 


A/ bh 


Mass of black hole 


10 


Yes 


9 ohs (degree) 


Disk inclination angle 


30, 60" 


Yes 


D (kpc) 


Distance 


10 


Yes 


rflag 


Self-radiation 







lflag 


Limb-darkening 







/col 


Spectral hardening factor 


1.7 


Yes 


K 


Normalization 


1 


Yes 


M 


Accretion rate 




No 


a 


Spin 




No 



b The values used in the fittings are the same as the value used in producing 
the spectra 

We show how the different SEDs in different models affect 
the measurement of the black hole spin. To achieve this, we 
create "fake" data sets from our spectra by convolving with a 
respo nse matrix of XMM-NEWTON using XSPEC (lArnaudl 
1996) version 12.5. The "faked" data sets are then fitted with 
the KERRBB model in XSPEC. 

The KERRBB model is a publicly available XSPEC model 
that can be used to calculate the multicolor disk (MCD) 
spectra of a relativistic, thin accretion disk around a Kerr 
black hole. The model settings in our fittings are summa- 
rized in Tab l e [D Our fitting strategy is similar to that of 
Shafe eetafl (|2006). The differences are that we turn off the 
limb-darkening and self-radiation, and we always set the spec- 
tral hardening factor f co \ to be 1.7. This is because we require 
the fittings consistent with the assumptions in our model. 
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Fig. 6. — The emergent spectra. The upper panels show the spectra when the disk is viewed at 30°, while the lower panels show the spectra 
at 60°. The left panels show the spectra from the non-spinning black holes (a = 0), the right panels show the spectra from the spinning black 
holes (a = 0.98). The black dashed lines correspond to the spectra from the relativistic standard disks (case A), the blue dotted lines correspond 
to the spectra form the relativistic slim accretion disks with no thickness (case B), and the red solid lines correspond to the spectra from the 
relativistic slim accretion disks with the inclusion of thickness effect (case C). The accretion rate takes the values: 1 /4M c( jd, 1 /2M c( jd, lM e dd> 
2M c( jd(left panels, a = 0) and l/16M c dd, 1/4A/ C dd, l/2A4edd> lM c dd, 2M e dd(right panels, a = 0.98) respectively. The black hole mass is 
taken to be 10M©. 
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The fittings are carried out by making our calculated spectra 
into a table model, then the faked data can be easily generated 
using the fakeit command. The faked spectra are then fitted 
to model KERRBB with the procedure described above. The 
fittings are robust in the sense that the different initial values 
of black hole spin a and the accretion rate M lead to almost 
the same results. 

Figure [9] shows the measured spin. For the moderate accre- 
tion rate, the effect of the heat advection alone leads to the 
overestimation of black hole spin, while for the high-enough 
accretion rate, this usually leads to the underestimation of the 
black hole spin. This is true for both the Schwarzschild and 
Kerr cases, and for the Kerr case these phenomena occur at 
a relatively lower accretion rate. The effect of the disk thick- 
ness will always lead to the under-estimation of black hole 
spin and is very significant when the accretion rate, the incli- 
nation or the spin, is large. Combining these two effects, in 
most cases, we generally underestimate the black hole spin. 
Only when the physical parameters are appropriate, we can 
obtain a overestimation of the black hole spin. 

7. DISCUSSIONS 

7.1. Shadowing Effect Enhanced by the Strong Gravity 

It is of theoretical interest to investigate what influences the 
effect of the self-shadowing. Without the light-bending effect, 
the self-shadowing effect is dependent on the maximum of 
H/r, which is mainly depend ent on the black hole accretion 
rate dMcClintock et al.1 120061) . However, it remains unclear 
how the shadowing effect influenced by the strong gravity. 

To gain physical insight into this question, we choose two 
systems with exactly the same luminosity but the different 
spin (one with M = 1 and a = 0.974, another with M=4 
and a = 0), and calculated the spectra. For both systems, we 
calculate the spectra both with and without the effect of light- 
bending. In the latter case, for self-consistency, the redshift is 
evaluated as 



Al/2 7 



1 



1/2 



1 ± 



7,5 



A 1 / 2 



(29) 



After that, to quantify the effect of disk self-shadowing, we 
compare the SED calculated with disk thickness and the SED 
calculated without disk thickness by defining the spectral dis- 
tortion as function of frequency as 

D. , . . / \ SED no t shadowed SED s h a dowcd 
istortion(j/j = — — , (30) 

^ED no t shadowed "T ^ED g h ac [owed 

where the SED not shadowed denotes the SED calculated by as- 
suming slim disk with no thickness, and SED s hadowed denotes 
the SED calculated including the disk thickness. 

The spectral distortion for different spin and different as- 
sumptions about light trajectories is plotted in Figure [10] For 
both the Schwarzschild and the Kerr black hole, including the 
light bending effect brings about significant distortion to the 
emergent spectra, and this distortion generally appears as at- 
tenuation of the flux at high energy part. This kind of distor- 
tion is more significant for Kerr black hole, because the light 
trajectories are more severely influenced by the strong gravity. 

7.2. Diversity in Spectral Evolution 

It was known that the spectral evolution of the X-ray bi- 
nary deviates from the trend predicted in the SSD model 



at both the low and high luminosities dDone et al.l 120071 : 
iKubota & Makishimal 120041) . Although the deviation at the 
low luminosity states is often inte rpreted as the evaporation 
dMeyer & Meyer-Hofmeisterll9 94) of the optically thick disk 
into the op tically thin ADAF (Advection Dominated Accre- 
tion Flow, iNaravan et all \l99% iNaravan & Yil [79951 [19911) . 
the deviation at the high luminosity states is still mysterious. 

The standard accretion disk model predicts the monotonic 
evolution of the disk peak energy with the luminosity, i.e., 
the higher peak energy corresponds to the higher luminosity. 
This kind of spectral evolution is consistent with the spectra 
evolution of X-ray binary at the medium luminosity states . 
However, beyond some critical luminosit y L cr j t the peak o f 
the disk SED shifts to the low energy band dDone et alj 2007). 
with the increase of the disk luminosity. This kind of spectral 
evolution is difficult to be understood in the SSD model. 

To tackle this problem, various attempts have been made, 
which include a change o f the color correction with the lumi- 
nosity (Da vis et alj20 06) as well as the effects of the optically 
thick advection. One common difficulty with these explana- 
tions is that they all predicts a single critical luminosity L cl -it, 
which is in contradiction with the diversit y of the critical lu - 
minosity observed in the different sources dDone et alj|2007l) . 

One natural explanation from our calculations is that the ab- 
normal of the spectral evolution is due to the self-shadowing 
of the accretion disks. The trend of the observed spectral evo- 
lution is consistent with our results as shown in Figure[6] and 
the diversity of the critical luminosity Ldisk/^edd in the dif- 
ferent sources might be due to the different inclination angle 
and black hole spin. According to our results, the critical lu- 
minosity becomes smaller for the larger inclination angle or 
the black hole spin. 

Our results also provide a diagnostic of the location of the 
shadowing medium. If the shadowing happens at the inner 
region of the accretion disks, the critical luminosity is depen- 
dent on black hole spin as well as the inclination angle. If the 
disk is shadowed by t he matter distant from the black hole, 
(e.g., equatorial wind; dDone & Da vis 2008)), the critical lu- 
minosity is dependent only on the black hole spin. 

7.3. The Spin of GRS 1915+105 

The spin of GRS 1915+105 is debated in the lit- 
erature. By f itting the disk blackbody component, 
McClintoc k et al.l (120061) foun d that the spin is near extreme, 
while Middleto n et al.l d2006l) found the spin is moderate, 
with a ~ 0.8. One major diffe rence between these two 
works is that the spin obtained by iMcClintock et alJ (120061) 
is based on the spectra fitting for t he relatively low luminos- 
ity (L < 0.3L odd ) states, while iMiddleton et all d2006l) set 
no strict restriction on t he luminosity states. Interestingly, 
IMcClintock et al.l (|2006) found that the measured spin be- 
comes smaller for the higher luminosity states. 

To show how the spin can be recovered by the fitting pro- 
cedure, we simulate the artificial data in ou r model which i s 
added with componization model comptt (Titarchuk 19941) . 
iron absorption edge, iron absorption line as well as galactic 
absorption, and then we fit the data with the same procedure 
as used in McClinto ck et alJ (l2006h. All the parameters are 
the same as those in McClintock et al. (2006). The only dif- 
ference in the artificial spectral fitting is that we replace the 
KERRBB2 model with KERRBB model, and set the spectral 
hardening factor to be 1.7. 

Our results are shown in Figure [TT] Without introducing 
the additional parameters, we are able to roughly reproduce 
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Fig. 7. — The angular dependence of the observed X-ray luminosity. As in Figure[6] the (black) dashed lines stand for the relativistic standard 
accretion disk, the (blue) dotted line stands for the dynamical thin accretion disk, and the (red) solid lines stand for the dynamical thick accretion 
disk. Different panels correspond to different black hole spin and accretion rate. The black hole mass is taken to be IOMq . 
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Fig. 8. — The angular dependence of the observed total luminosity. The symbol of the lines are the same as in Figure|7] 
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FIG. 9. — Measured spin for different models. The upper panel shows the 
result for a non-spinning black hole with a = 0, and the lower panel shows 
the result for a spinning black hole with a = 0.98. The dashed lines corre- 
spond to the thin Keplerian case (case A), the dotted lines correspond to the 
slim accretion disks with no thickness (case B), and the solid lines correspond 
to slim accretion disks with thickness effect (case C). Different line thickness 
represent different viewing angles: Thin lines for 8 = 30°, and thick lines 
for 6 = 60°. 
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FIG. 10. — Spectral distortion (defined by Equation l30t as a function of the 
photon energy. The solid lines standard for the results with the effect of the 
light bending, while the dotted lines for the results without the light-bending 
effect. The red lines stand for the Kerr black hole with a = 0.98, and the 
black lines for the black hole with a = 0. In all the cases, the disk luminosity 
is the same. 



the trend found in iMcClintock et alj (120061) . At both the low 
(around 1/4 of Eddington luminosity) and high (close to 1 
Eddington luminosity) luminosity states, our simul a ted re - 
sult agrees well with the result of McClint ock et al.1 (2006). 
When the disk luminosity is low enough, according to the 
spectral model we have developed, the spin determination by 
IMcClintock et al.l d200 6) is reliable. When the luminosity is 
high enough, the disk shadowing makes the measured spin 
lower than the actual value. This luminosity dependence of 
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FIG. 1 1 . — Measured spin of GRS 1915+105 as a function of the disk lu- 
minosity. The green line and dots are the results calculated in our model, 
including the effects of the heat advection and the disk s elf-shadowing. The 
red dots are the spin obtained in McClintock et al. 12006) by fitting the ob- 
served SED. 

the measured spin is shown up both in our calculation and in 
the observational data. When the disk luminosity is medium 
(around 10~ 02 Eddington luminosity), our results begin to 
diverge from the data. This may due to the photon scat- 
tering by the wind fr om the accretion disk as suggested by 
iDone & Davis! d2008|). However, as our s pectra fitting is based 
on the model of McClinto ck et all (120061) . so it is still possible 
that different spectral model (e.g. a low tem perature comp- 
tonization component, Middl eton et al.ll2006T) may lead to a 
different result. 

8. CONCLUSIONS 

In this work, we study the structure of the relativistic slim 
accretion disks and their emergent spectra. Besides the effect 
of the photon trapping, we find that when the accretion rate 
is appropriate, the energy trapped in the outer regions of the 
accretion disks can be re-radiated in the inner regions of the 
accretion disks. Another main result of this work is that we 
check the validity of the relativistic SSD. We find that at low 
accretion rate, such as M ~ 1/4 for a — and M ~ 1/16 
for a = 0.98, the relativistic SSD is still valid to describe the 
disk structure, and the effect of heat advection will not signif- 
icantly influence the global structure of the disk. At relatively 
high accretion rate such as M = 2Af c dd, these effects will 
significantly distort the emergent spectrum, and hence affect 
the measurement of the black hole spin. We also find that 
when the accretion rate is high, the effect of the disk self- 
shadowing plays an even more important role than that of the 
heat advection alone. 

The emergent luminosity from the accretion disks is angu- 
lar dependent, and the degree of anisotropy is dependent on 
the accretion rate. In this work, we investigate the angular 
dependence of the observed luminosity. Due to the physi- 
cal similarities between the accretion disks of different black 
hole mass, our result can be applied to estimate the radiation 
anisotropy of AGNs. 

To check in principle the effect differences of spectra on es- 
timation of spin, we create "fake" data sets with the simulated 
data, and fit them with the KERRBB model using XSPEC. At 
the relatively low accretion rate, our fitting results suggest that 
the simple Keplerian model is reliable. At the relatively high 
accretion rate, our fitting results suggest that self-shadowing 
effect will lead to significant under-estimation of the black 
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hole spin. By making our calculation relevant to the case of 
GRS 1915+105, without introducing additional parameters, 
we are able to reproduce the decrease of m easured spin with 
the inc reasing luminosity, which is found in McCl intock et alJ 
(2006). This finding suggests that disk self-shadowing signif- 
icantly shapes the spectra of GRS 1915+105. 

There are several effects which we do not include in our 
calculations. It was claimed that the disk may have no time t o 
relax to the thermodynamic equilibrium (Ohsuga et al. 2002), 
this effect will make our estimation of emergent flux (Equa- 
tion ( [151 )) invalid and make the disk thickness estimation in- 
accurate. In our calculations of the disk SED, we use a sim- 
plified model of the local spectrum and neglect the effect of 
limb-darkening and returning irradiation. Assuming a con- 
stant color correction may make our calcula ted spectra in- 
adequ ate to be compared to the observation dDone & Davis! 
2008), so we focus on analyzing the physical effects we have 
introduced. The limb-darkening effect will influence the spec- 
tra and the angular dependence of the emergent luminosity, 



while the returning irradiation will have some effects on both 
the disk structure and the emergent spectra. As lLi et ail d2005l) 
have calculated, these two effects will contribute several per- 
cents to the final results. Despite of these shortcomings, our 
results and the trend of the spectral evolution are still robust 
when the inclination angle is not too large. 
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